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We present self-consistent calculations of the spin-averaged pair distribution function g(r) for 
a homogeneous electron gas in the paramagnetic state in both three and two dimensions, based 
on an extension of a model that was originally proposed by A. W. Overhauser [Can. J. Phys. 
73, 683 (1995)] and further evaluated by P. Gori-Giorgi and J. P. Perdew [Phys. Rev. B 64, 
155102 (2001)]. The model involves the solution of a two-electron scattering problem via an effective 
Coulombic potential, that we determine within a self-consistent Hartree approximation. We find 
numerical results for g(r) that are in excellent agreement with Quantum Monte Carlo data at low 
lO ■ and intermediate coupling strength r s , extending up to r s m 10 in dimensionality D = 3. However, 

the Hartree approximation does not properly account for the emergence of a first-neighbor peak at 
stronger coupling, such as at r s = 5 in D — 2, and has limited accuracy in regard to the spin-resolved 
components <7n( r ) an d 5n( r )- We also report calculations of the electron-electron s-wave scattering 
| length, to test an analytical expression proposed by Overhauser in D = 3 and to present new results 

S-H , in D — 2 at moderate coupling strength. Finally, we indicate how this approach can be extended 

^ ' to evaluate the pair distribution functions in inhomogeneous electron systems and hence to obtain 

improved exchange-correlation energy functionals. 
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PACS numbers: 71.10.Ca - Electron gas, Fermi gas 



O \ 71.45.Gm - Exchange, correlation, dielectric and magnetic functions, plasmons 



I. INTRODUCTION 

Many of the electron-electron interaction effects in simple metals and semiconductors can be understood by reference 
to the homogeneous electron-gas model. A central role is played by the electron-pair distribution function g(r), which 
describes how short-range exchange and correlations enter to determine the probability of finding two electrons at a 
relative distance r.0 The mean potential energy of the electron gas can be calculated from g(r) and hence, through 



an integration over its dependence on the coupling strength, the total energy as well. Thus, an accurate knowledge 
of this function is crucial for applications of density functional theory! (DFT) in various schemes that have been 
I . proposed to transcend the local density approximation (LDA) in the construction of exchange and correlation energy 
functionals .! 

A precise definition of g(r) is through the average number of electrons contained in a spherical shell of radius 
' r and thickness dr centered on an electron at the origin, which is given by ng(r)ilD rD1 dr where D is the space 
dimensionality, n the electron density, and fi^ the solid angle in D dimensions (with f2 2 = 2-7T and f2 3 = 47r). 
In fact, g(r) is the average of the distribution functions for parallel- and antiparallel-spin electron pairs, g(r) = 
[i?Tt( r ) + fftl( r )]/2- While this spin-average reflects the charge-charge correlations in the electron gas, the spin-spin 
correlations are described by the difference distribution function gd(r) — [.9u( r ) — ffTl( r )]/2- 

Early calculations of the pair distribution functions in the 3D electron gas were based on the use of a Bijl-Jastrow 
correlated wave function for the ground state! and on exploiting the fluctuation-dissipation theorem for a self-consistent 
determination from the charge and spin response functions. □ These early results were validated with the advent of 
the Quantum Monte Carlo (QMC) techniques,0 which have produced a wealth of accurate data on correlation and 
response functions over a wide range of coupling strength for both the 3Z)B! and the 2fM case. The QMC data have 
in turn stimulated a number of further theoretical studies. 
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Here we are specifically concerned with the approach proposed by OverhauserS for the evaluation of the value 
g(0) of the pair distribution function at contact. While g-f-f(O) vanishes on account of the Pauli principle, <7fj,(0) 
is determined by the two-body scattering events and in Overhauser's model was obtained from the solution of an 
effective Schrodinger equation for the relative motion of two electrons with antiparallel spins. Earlier work had 
demonstrated an exact cusp conditionQ relating the logarithmic derivative of Qj | (r) at contact to the Bohr radius, 
and had emphasized the importance of the electron-electron ladder diagramst3 in evaluating g^i(0). The approach 
of Overhauser has subsequently been used to evaluate <?|j,(0) in the 2D electron gas,El and has been extended by 
Gori-Giorgi and PerdewEj to evaluate g(r) at finite r in 3D through an accurate numerical solution of Overhauser's 
two-body Schrodinger equation. Their results are in good agreement with QMC data in the short-range part of g(r). 

In the present work we further develop this approach to the evaluation of the pair distribution functions by (i) adopt- 
ing a self-consistent Hartree scheme for the determination of the effective potential entering the two-body Schrodinger 
equation; and (ii) carrying out calculations for both a 3D and a 2D electron gas. The paper is organized as follows. 
Section [n] presents the essential theoretical background and Section [II reports and discusses our numerical results. 
In Section |w] we indicate how this approach could be extended to evaluate pair correlations in an inhomogeneous 
electron gas and hence to obtain improved exchange-correlation energy functionals for applications of DFT. A brief 



summary concludes the paper in Section IV 



II. ESSENTIAL THEORY 



Following the work of OverhauserEEl and of Gori-Giorgi and PerdewEfl, we aim to solve the two-electron scattering 
problem in some effective interaction potential V(r) (assumed to be independent of the spin state of the electron 
pair) in order to determine the wave functions f[' K (r, R) and ^ R) for the parallel and antiparallel spin states. 
Here, r = ri — r2 and R = (ri + r2)/2 are the relative coordinate and the center-of-mass coordinate of the pair, the 
conjugate momenta being k and K. The spin-resolved pair distribution functions g aa > (r) can then be obtained as 



wMH(l*Mc(r,R)l 2 M pW 



(1) 



where ((■■■)e)p(k) represents averages over the scattering angle and over the probability p(k) of finding two electrons 
with relative momentum k in the electron gas, 



(2) 



((G)e) P (k) = (l/nu) / dk P (k) / dn D o 



The function p(k) can be calculated from the momentum distribution 



using 
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The prefactor in front of the integral is obtained from the condition that p(k) should integrate to unity. 
The pair wave functions can be written in the form of angular- momentum expansions, 
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where Af\6) = (21 + l)P^(cos6») in terms of the Legendre polynomial P e (x) and Af\ff) = 2cos(£9). In Eq. (g) 
the sum over £ runs over odd integers, because the spinor associated with the |T state is symmetric and hence the 
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coordinate part of the wave function has to be antisymmetric under exchange. The spinor associated with the fj. state 
has instead no definite symmetry and hence the sum in Eq. (||) runs over all integer values of £. Upon performing 
the angular average in Eq. (0) one obtains 
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1 and B^p = 2. Finally, the Schrodinger equation for the wave function (r) is 
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+ F(r)U^(r) = 
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where /i = m/2 is the reduced mass of the electron pair, L| = hf £{£ + 1) with i = 0, 1, 2, 



and L| = ?i 2 ^ 2 with 



0, ±1,±2, ... (the negative values of £ are accounted for by the choice B^ = 2 in Eqs. (||) and ([[])). In solving 



Eq. (|S|) we impose as a boundary condition that ^/^{r) tends asymptotically to the free solution (i.e. the one which 
is obtained by setting V(r) = 0) except for a phase shift. 

The form of Eq. (||) ensures that the relation g-\-\{§) = is satisfied, since all functions ^/^(r) vanish at the origin 
for I ^ 0. The cusp condition on g\i[r) reads d\nn^ ^(r) / dr\ r= o — 1/as in D = 3 and din g-\i(r) / dr\ r= o = 2/as in 
D = 2: following the argument given by KimballM, it will be satisfied if V(r) tends to the bare Coulomb potential 
for r — > 0. It is also easily seenli-3 that the Hartree-Fock results for g^(r) and S'-fj.(r) are recovered if V(r) is set to 
zero in Eq. (0). Given a general scattering potential V(r), there is no guarantee that the charge neutrality condition, 



d D r [g{r) - 1] 
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(9) 



is satisfied. We have numerical evidence that Eq. (||) is fulfilled for our choice of V(r) (see subsection A). The accuracy 
with which this happens depends in practice on the numerical solution of Eq. (||) and thus on the number of angular 
momentum states that are included in Eqs. (|) and (1)0 

Before proceeding to present our choice for the potential V(r) in Eq. (||), we report the expressions for the probability 
function p(k) introduced in Eq. (^|). Using the momentum distribution n(k) — 9(kp ~ k) for the free Fermi gas, one 
obtains 
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in D = 30 and 
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m D = 2.0. We have also evaluated p(k) for interacting electrons using QMC data on the momentum distribution 
from the work of Ortiz and Ballonei in 3D and of ContiEll in 2D. However, this led to only small changes in the pair 
distribution functions reported in Sect. Ill with the help of Eqs. (|Io|) and (111]). 



A. Hartree potential and self-consistency 

Overhauser's proposaS for calculating g(0) in the 3-D case was to approximate V(r) by the electrical potential of 
a model consisting of an electron at the origin and a neutralizing sphere of uniformly distributed charge with radius 
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r s a,B = (47rn/3) -1 / 3 . Gauss's law ensures that V(r) vanishes outside the sphere, and an approximate solution of 
the scattering problem could be obtained by an iterative procedure. In the work of Gori-Giorgi and PerdewEl the 
same model for V(r) was adopted to evaluate g(r), but a full solution of the Schrodinger equation (||) was achieved. 
In the 2D case with e 2 /r interaction, on the other hand, the potential outside a uniformly charged disk of radius 
r s a,B = (m) -1 ' 2 with an electron at its center does not vanish, since the electrical force field extends outside the 
plane in which the electrons are moving. A more rehned model is therefore necessary.El 

Here we approximate V(r) in Eq. (||) by the Hartree potential due to the whole distribution of electrical charge and 
evaluate it with the help of Poisson's equation. More precisely, in the 31? case V(r) is taken to satisfy the equation 

A r V(r) = -4vre 2 [(5(r) + n(g(r) - 1)] , (12) 

where A r is the radial Laplace operator. The appropriate Poisson equation for the Hartree potential Vn(r, z) in the 
2D case, with r the radial distance in the electron plane and z the vertical distance from the plane, is 

(A r + ^)Vk(r, z) = -Ane 2 [5{r) + n(g(r) - l)]6(z) (13) 

and what is needed is V(r) — Vn(r, 0). Evidently, by solving Eq. (||) in conjunction with Eq. (12) or Eq. ( |l3| ) we 
obtain a self-consistent determination of the effective potential and of the radial distribution function. 

The solution of Eq. ( [12] ) and Eq. ([l3]) is easily obtained by introducing Fourier transforms. We define the structure 
factor S(k) through the relation 

S(k) = 1 + n J d D r [g(r) - 1] exp (-ik ■ r) . (14) 

It is then easily seen that the Fourier transform of V(r) (V(k), say) is given by 

V{k) = v{k) S(k) (15) 

where v[k) — 4ire 2 /k 2 in D — 3 and v(k) — 27re 2 /k in D = 2. To prove Eq. ( [l5|) in the 2D case, we notice that the 
Fourier transform of Vn(r, z) from Eq. ( |l3| ) is 

~ 4-7TP 2 

v ^ k > k ^ = WTk~l s{k) (16) 

so that 

y H (fc,z)^ J^ r e^ z V K (k 1 k z ) = ^fe- k ^S(k), (17) 

and V(k) = V R (k,z = 0) = 2ne 2 S{k)/k. 

Let us examine the asymptotic behaviours of the effective potential V{k) given in Eq. (|l5|). At large momenta the 
structure factor S(k) tends to unity, so that V(k) — > v(k) and V(r) tends to the bare Coulomb potential in the limit 
r — > 0. As already remarked, this property ensures that the cusp condition on g(r) is satisfied. In the limit k — > 0, on 
the other hand, charge neutrality ensures the validity of the plasmon sum rule, which may be written in the form 

lim5(fc) =e k /fkJ pi (18) 

where Ek — h 2 k 2 /2m and lu p \ is the leading term in the plasmon dispersion relation, given by u> p \ = (Anne 2 /m) 1 / 2 in 
3-D and by w p i = (2ime 2 k/m) 1 / 2 in 2D. Therefore, 

lim V"(jfe) = hu p i/2n (19) 

in both dimensionalities. That is, the Fourier transform of our choice for the effective scattering potential tends in 
the long-wavelength limit to a constant in D = 3 and to zero with a k 1 / 2 law in D = 2. 
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B. s-wave scattering length 



We complete this discussion by showing how the s-wave scattering length a sc (r s ) can be evaluated from the numerical 
solution of the electron-electron scattering problem. The s-wave phase shift 8o(k) is introduced through the large- 
distance behavior of the two-particle scattering state with £ = at fixed momentum fc, 

<kM ~ cos [kr - (D - if- + 6 (k)} . (20) 

From So(k) the scattering length is obtained by the requirement that the wave function outside the range of the 
potential should vanish at r ~ a sc . In D = 3 this yields from Eq. ( p(i|) the well-known relation 

a sc (r s ) = lim (21) 

A simple analytical expression for a sc (r s ) in the 3-D electron gas is available from the work of Overhauser0, 

/ s rJlO 
a S c{r s ) = r s a B — — — • (22) 
1 + 3r s /8 

The introduction of the concepts of effective range and scattering length is much more delicate in D = 2. In the 
work of Verhaar et aZ.EH the H-H atomic scattering problem was analyzed in detail. The most appropriate form of the 
outer wave function in 2D is 

*o*( r ) ~ exp [»*„(*)] cos [S (k)]{3 Q (kr) - tan [S (k)] N (fcr)} , (23) 

in terms of the Bessel functions Jq(x) and Nq(x). The relation between scattering length and phase shift thus is 
tan[<5o(fc)] = Jo(fca sc )/No(fca sc ), taking at low energy a form which is the same as for the scattering of two hard 
spheres of radius a sc , 

cot [6 {k)) = (2/vr) [ 7 + In (ka sc /2)} + o(fc 2 ) (24) 



where 7 is the Euler's constant, 7 = 0.577215665 • • • . In Sect. Ill we determine the scattering length for the 2D 
electron gas by fitting the expression in Eq. (^) to the phase shift obtained from the asymptotic behavior of the 
two-electron wave function as a function of k at low energy. An analytical expression for a sc in the strong-coupling 
limit has been given by Polini et al^. 



III. NUMERICAL RESULTS 

We report in Figures |l| to ^| our numerical results for the pair distribution functions and for the self-consistent 
scattering potential. Starting with the 3D system, Figure |l| shows the spin-averaged g(r) at r s = 10, in comparison 
with the QMC data reported by Ortiz et alM The results obtained by Gori-Giorgi and PerdewEl within the same 
theoretical scheme, but with Overhauser's choice for the effective potential V(r), are also shown in Figure [l]. It is 
seen from the Figure that both choices of V(r) yield excellent agreement with the QMC data for the short-range 
part of the electron-electron correlations, up to r/(r s as) w 1. However, the self-consistent calculation based on the 
use of the Hartree potential becomes definitely superior at large distance, where (at this intermediate value of the 
coupling strength) it continues to be in excellent agreement with the data. The Hartree potential at self-consistency 
is shown in Figure |^, for both the 3D and the 2D case. We may also remark that the cusp condition is satisfied by 
our numerical results in both cases. 

With further increase in the coupling strength the pair distribution function from the QMC work starts developing 
a first-neighbor peak, and this behavior is not reproduced quantitatively by the theory. This is illustrated in Figure 
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H for the electron gas in D = 2. As is well known, the reduction in dimensionality enhances the role of the electron- 
electron correlations: this is also clear from comparing the scattering potentials in the two panels in Figure @. While 
in 2D the self-consistent theory remains quite accurate at moderate coupling strength as is shown by the comparison 
with the QMC data of Tanatar and Ceperleyi at r s = 1 in the left-hand panel in Figure |3|, quantitative differences 
from the QMC data are seen to arise at r s — 5 (right-hand panel in Figure ||) . 

The other aspect of the theory that needs testing concerns the quality of its predictions in regard to the spin- 
resolved pair distribution functions. This point is examined in Figure |J for the 31? system at r s — 5 and 10, using 
the QMC data of Ortiz et al. as analyzed by Gori-Giorgi et aL.i The discrepancies between theory and simulation 
are reasonably small at these values of the coupling strength. It is evident from the Figure that these discrepancies 
largely cancel out in taking the spin avarage, but will be magnified when one calculates the difference distribution 
function gd(f)- Similar theoretical results are shown in Figure]^ for the 2D system. 

Finally, in Figures |and@ we report our results for the s-wave scattering length as a function of coupling strength 
r s . For the 3D system the simple analytical formula obtained by OverhausercJ and reported in Eq. ( p2| ) is seen in 
the left panel in Figure |^ to give a very good account of our results. Results for the 2D electron gas at moderate 
coupling strength are shown in the right-hand panel in Figure while Figure [?] shows how they have been obtained 
by fitting the expression ( p4| ) to our numerical results for the phase shift in s-wave as a function of momentum at low 
momenta. It may be remarked that the magnitude of So(k) in the present electron-electron scattering in 2D is smaller 
than that in the H-H scattering problem studied by Verhaar et al&3 by a factor of about 2. This yields, however, huge 
differences in the magnitude of the scattering length. Our results in the right-hand panel in Figure ^ should thus be 
regarded as being very sensitive to the details of the theory and hence as having limited quantitative value. 

IV. EXTENSION TO INHOMOGENEOUS SYSTEMS 

In this section we briefly indicate how the approach that we have presented in Sect. || could be extended to deal with 
the pair distribution function in an inhomogeneous electron system, subject to an external scalar potential I4 x t(r). 

The electron-electron correlations are described in such a system by an inhomogeneous pair distribution function, 
g(r, r') say. The exchange and correlation energy functional is given by 

E xc = \ J d D r J d D r> »(r)n(r') [g(r, r') - 1] »(|r - r'|) , (25) 

where n(r) is the inhomogeneous electron density and g(r,r') is obtained from ofr,r') by an integration over the 
coupling strength at fixed n(r) (see for instance the book by Dreizler and Grossa). The calculation of g(r, r') by 
means of a two-electron scattering problem remains in this case a genuine two-body problem. It requires for each spin 
state solution of the equation 



h 2 

-— (A ri + A r2 ) + y cxt (ri,r 2 ) + V(r u r 2 ) 
Am 



$e 12 (ri,r 2 ) = ei 2 $ £l2 (r 1; r 2 ) , (26) 



where A r is the Laplace operator, V^ xt (ri,r 2 ) = K xt (ri) + V^ xt (r 2 ), and F(i"i,r 2 ) is the effective electron-electron 
potential. The Hartree approximation on the effective potential, leading to 

V(r 1 ,r 2 ) = i;(|r 1 -r a |)+ J d D r' [n(r') g(r u r') - n b (v')\ v{\v' - r 2 |) (27) 

provides an approximate self-consistent closure of the problem transcending the usual LDA or other approaches that 
appeal to the exchange-correlation hole of the homogeneous electron gas. In the previous equation (f27|), nfe(r') is the 
density of the background. 

Finally, a relation between g(r, r') and the two-electron scattering states $ei 2 ( r : r ') is needed. This relation can in 
general be written in the following form: 



6 



<w(r,r')= Yl r^'|$ ei2 (r,r')| 2 , 

ei2 , occup. 



(28) 



where the sum runs over all the occupied levels ei2- The appropriate degeneracy factor T™ 2 for the eigenvalue e 12 is 
zero in the case a — a' , if $ei2( r i r ') i s symmetric under the exchange r <-> r'. 

For building the exchange and correlation energy functional in Eq. (|2^) , one needs to calculate the pair distribution 
function at each given coupling strength A by repeating the procedure outlined above with v(\ri — i"2 1) = A/|ri — r2 1 . 
The density profile n(r) is needed at full coupling strength and requires to be obtained by a parallel DFT procedure. 
The pair distribution function at full strength is, in itself, a very interesting quantity. 

As an example of application of such scheme, we would like to mention the problem of a finite number of electrons 
confined in a quantum dotS. In this case the eigenvalues £12 are discrete and the summation procedure reported in 
Eq. ( p8| ) corresponds to filling the lowest energy states with all the available electrons. 

Another problem of interest is represented by a system of electrons confined in a quantum well, in cases where 
the electron dynamics is important also in the growth direction so that the confinement cannot be handled by a 
simple reduction to 2D. Compared to the previous example, the summation procedure is in this case more involved. 
The difficulty comes from the fact that the motion in the transverse direction is free, and this implies that for 
each subband in the growth direction there is a dispersion in the transverse direction associated with the in-plane 
momentum. Examples of such physical systems, which are in a sense intermediate between 3D and 2D, have been 
discussed for instance by Ullrich and Vig nalei and by Luin et alB. 

More generally, the possibility of testing theories of the exchange-correlation hole and of the exchange-correlation 
energy density in strongly inhomogenous electronic systems is opening up through novel applications of variational 
QMC methods!!. 



V. SUMMARY 

In this work we have proposed an extension of Overhauser's model for the electron-electron correlations in the 3-D 
electron gas on the basis of a self-consistent Hartree approximation for the electron-electron scattering potential. We 
have confirmed that the model is quite accurately describing the short-range part of the exchange-correlation hole, as 
already demonstrated by Gori-Giorgi and PerdewB, and shown that it can be usefully extended to cover the full range 



of interelectronic separation over an appreciable range of values of the coupling strength. As already noted in Sect. Ill, 
the accuracy of the present Hartree approximation is limited to the range of values of r s below the development of 
a first-neighbor peak in the pair distribution function. We have also shown that the original proposal of Overhauser 
yields a very accurate analytical formula for the electron-electron scattering length, and indicated how this approach 
could be extended to deal with the exchange-correlation hole in an inhomogeneous electron gas, leading perhaps to 
more accurate descriptions of the exchange-correlation energy functional. 

We have examined the usefulness of this approach in describing the exchange-correlation hole in the 2D electron 
gas and the spin-resolved pair distribution functions in the 31? electron gas. As may be expected, the Hartree approx- 
imation is quantitatively useful in 2D over a more limited range of coupling strength and has more limited accuracy 
in regard to the splitting of the exchange-correlation hole into its parallel- and antiparallel-spin components. With 
regard to these spin-resolved pair functions, however, we feel that it is quite remarkable that a Hartree approximation 
should already work as well as it does in our calculations. We may hope that major changes will not be needed to 
explicitly include many-body exchange corrections in the electron-electron scattering potential. 

We have also given some attention to the determination of the electron-electron scattering length in the 2D electron 
gas. As is generally true for 2D systems, the logarithmic dependence of the s-wave phase shift on ka sc at very low 
energy makes this problem a very delicate one. We have shown that the proposal made by Verhaar et aZ.EH in the 
context of an atom-atom scattering process is just as useful in regard to electron-electron scattering. 
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We are currently examining how the Hartree scattering potential may be supplemented by explicit inclusion of 
corrections accounting for many-body exchange and correlations, with the main aims of (i) improving the quantitative 
account of the spin-resolved components of the exchange-correlation hole, and (ii) studying how a first-neighbor shell 
emerges with increasing coupling strength. 
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FIG. 1. The pair distribution function in a 3D electron gas at r 3 — 10, as a function of r/(r s a»i. The results of the 
self-consistent Hartree approximation (fulli4ine) are compared with QMC data (crosses, from Ortiz et al.u) and with the results 
of calculations by Gori-Giorgi and Perdewi-1 (dotted line). 
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FIG. 2. Self-consistent scattering potential V(r) in the Hartrae-iapproximation (in units of h 2 kp/m), as a function of jp//r s as) 
in 3D (left) and in 2D (right). The Overhauser potential in 3Dcj at r s = 10 and the result of the work of Polini et al&3 in 2D 
at r s = 5 are reported for comparison (dotted lines). 
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FIG. 4. The spin-resolved pair distribution functions in a 3D electron gas at r a = 5 and r a = 10. The resists of the 
self-consistent Hartree approximation (full lines) are compared with the QMC data (crosses, from Gori-Giorgi et alH). 
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FIG. 5. The spin-resolved pair distribution functions in a 27? electron gas at r s = 1 and r s — 5, from the self-consistent 
Hartree approximation. 
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FIG. 6. The s-wave scattering length as a function of r s , in units of the Bohr radius. Left panel: the result of the ««lf-consistent 
Hartree approximation for the 3D electron gas (full line) is compared with the analytical result of OverhausertJ in Eq. jj^ ) 
(broken line). Right panel: results from the self-consistent Hartree approximation for the 2D electron gas. 
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FIG. 7. The s-wave phase shift in the 2D electron gas at various values of r s , as a function of reduced momentum k/kF at 
low momenta. Full lines: from the numerical solution of the self-consistent Hartree approximation; broken lines: from Eq. ( fj^ ) 
with the values of the scattering length shown in the right-hand panel in Figure 0. 
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